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ABSTRACT 

We report here a study of gas, dust and star formation rates (SFRs) in the molecular cloud complexes 
(MCCs) surrounding the giant HII region RCW 106 using 12 CO and 13 CO (1-0) data from the Three- 
mm Ultimate Mopra Milky way Survey (ThrUMMS) and archival data. We separate the emission 
in the Galactic Plane around l = 330°-335° and b = —1°-1° into two main MCCs: the RCW 106 
(Vlsr = —48kms _1 ) complex and the MCC331-90(Vlsr = — 90kms -1 ) complex. While RCW 106 
(M~ 5.9 x 10 6 M e ) is located in the Scutum-Centaurus arm at a distance of 3.6 kpc, MCC331-90 (M~ 

2.8 x 10 6 M©) is in the Norma arm at a distance of 5 kpc. Their molecular gas mass surface densities 
are ^ 220 and ^ 130M© pc -2 , respectively. For RCW 106 complex, using the 21cm continuum 
fluxes and dense clump counting, we obtain an immediate past (~-0.2 Myr) and an immediate future 
(^+0.2 Myr) SFRs of 0.25^0.023 M©, yr _1 and 0.12 db 0.1 M© yr -1 . This results in an immediate past 

SFR density of 9.5^9 M© yr _1 kpc -2 and an immediate future SFR density of 4.8lgJ M© yr -1 kpc -2 . 

As both SFRs in this cloud are higher than the ministarburst threshold, they must be undergoing 
a ministarburst event although burst peak has already passed. We conclude that this is one of the 
most active star forming complexes in the southern sky, ideal for further investigations of massive star 
formation and potentially shedding light on the physics of high-redshift starbursts. 

Subject headings: stars: formation, stars: protostars, ISM: clouds, ISM: structure, ISM: HII regions 


1. INTRODUCTION 

Studying the earliest phases of star formation involves 
examining the morphological and kinematic structure of 
molecular clouds and the transformation between differ¬ 
ent states of materials. Massive star formation often oc¬ 
curs in molecular cloud complexes (MCCs) with radius 


70 p_c_or more (e.g., W43; e.g., Nguyen Luong et al. 
2011a Motte et al.|2014 ), Cygnus X with radius ” 80 pc 
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(e.g., Schneider et al. 2006 Motte et al. 2007), and the 
Cent ral M olecula r Zone (CM Z) with radius ~ 180 pc 
(e.g., [Miyazaki & Tsuboi|2000|| Jones et al.|2012|). MCCs 
associated with massive star formation reside mainly 
along the mid-plane of the Galaxy, may be distant and 
heavily obscured in the visible regime. Though not be¬ 
ing obscured in the radio wavelength, Large-field radio 
surveys of MCCs generally yield data sets with a sub¬ 
stantial field of view but insufficient resolution to con¬ 
duct case-by-case studies of particular clouds. This is 
particularly true for low resolution surveys such as the 


carbon monoxide (CO) survey of th e Milky Way (Bronf¬ 


man et al. 


_B y 

Dame et alj 2001). While they catalog 


molecular clouds, the structure and properties of each 
cloud cannot be fully resolved. However, recent high- 
angular-resolution spectroscopic surveys of the Galactic 
plane such as the Galactic Ring Survey (GRS; Jackson 


et al. 2006| ) and the Thr ee-mm Ultimate M opra Milky 


way Purvey (ThrUMMS; |Barnes et al.||2013|) help to ad¬ 
dress this problem. 

Of particular interest here is the MCC of gas surround¬ 
ing the bright giant H n region RCW 106, which was dis¬ 
covered in t he Hg emission line s urvey of the southern 
Milky Way (Rodgers et al. 1960). The giant H n re¬ 
gion RCW 106 hosts a cluster with mass ~ 10 3 M© and 
a Lyman continuum photon emission of 10 50 s -1 , likely 
originating from dozens of O-type stars (M > 8 M©). 
It resides in the Scutum-Centaurus arm at a d istance 
of ~ 3.6 kpc from the Sun (Russeil et al.|2005]). This 
makes it an excellent laboratory for testing molecular 
cloud formation mechanisms in response to galactic dy¬ 
namical processes such as spiral density waves. 
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Table 1 

Observation parameters of the data used in this study 


Instrument 
or Survey 

Tracer 

Frequency 

(GHz) 

Ang. Res 

o 

Velocity 
(km s -1 ) 

Av res 
(km s -1 ) 

la rms 

CfA 

12 CO 

115.27 

510 

-165 to 165 

1.3 

0.1 K km s -1 

ThrUMMS 

12 CO 

115.27 

72 

-140 to 10 

0.33 

~ 3 K km s -1 

ThrUMMS 

13 CO 

110.20 

72 

-140 to 10 

0.34 

~ 1 K km s -1 

SGPS 

H I 

1.420 

132 

- 

- 

1.6 K km s -1 

SGPS 

21 cm cont. 

1.420 

132 

- 

- 

1 mjy beam -1 

MIPS /Spitzer 

24 ym 


6 

- 

- 

0.01 MJy sr -1 

PACS / Herschel 

160 ym 


12 

- 

- 

0.08 MJy sr -1 

SPIRE /Herschel 

500 ym 


37 

- 

- 

1.2 MJy sr -1 


Moreover, the parent RCW 106 MCC is also the site 
of ongoing star formation. It is centered on l ~ 333°, 
b 


Ulsr ^ —50 kms 1 (Bains et al. 

2006 

).Its substructures 

have been subject to a variety oj 

t spectral line and con- 


tinuum studies, progressively unravelling the star forma- 
tion and evol utiona ry history of the regio n. Specifically, 


(Lynga 1964) and (Urquhart et al. 2007) measured the 
impact of UV emitted by various nearby OB stars that 
ionize the surrounding environment. 1.2 mm continuum 
maps revealed 95 clumps with masses ranging from 40 to 
10 4 M^, some of which have infrared (IR) coun terparts 


suggesting embedded star formation (Mookerjea et al 
2004). Higher density tracer such as 13 CO, CS, HCO + , 


HCJNI, and HNC emission lines show prominent velocity 
features centered on the U LSR of RCW 106 that coincide 
with 1.2 mm d ust clumps and are also sites of massive 


star forma tion (Bains et al. 2006 Wong et al. 2008 Lo 


et al.|2009| ) . A recent N H 3 survey toward dense clumpsol 


the RCW 106 complex revealed a large sample of col d col¬ 
laps ing clumps that are potentially forming stars (Lowe 
etah|20l4| ). 


The two densest and most massive sites of o ngoing 


et al. 


pact H ii ( I 

2005, [2006) 


star formation core are MMS5 and MMS68 (Mookerjea 
|2004|). They are associated wit h the ultracoim 
4 ii (UCHu) regions G333.6-0.2 (|Fujiyoshi et al. 
and IRAS 16164-5096, respectively! High 
resolution adaptive optics observations with OZ'3 resolu¬ 
tion 1000 AU) combined with mid-IR spectroscopy 
indicate the presence of mid_to late-typ e O stars that 
might still be gaining mass ( Grave et al.||2014 ). 


gas and dust at cloud scales of 0.1-10 pc ( Mookerjea 

et al. 2004; Bains et al. 2006; Wong et al. 2008; |Lo et al. 

2009 

), dense clump scales of 0.1-1 pc (Lo et al. 

to 

o 

o 

2011 

), or embedded O stars subparsec scales oi 

t even 

smal 

ler structure (Kumar 2013 Grave et al. 2014). In 


50 — 100 pc), which is the largest scale of molecu¬ 
lar cloud in the Galaxy but the smallest scale that can 
be probed in nearby galaxies. On this scale, RCW 106 
is a massive filamentary molecular cloud complex with 
an aspect ratio larger than 2: see Figure [5j Our goal 
is to study the global properties of the molecular clouds 
surrounding RCW 106 and establish a relationship be¬ 
tween the large scale structure of these clouds and their 
enhanced star formation. 

The various data sets used in this paper are described 
in Section [2] We present the two star-forming MCCs 


identified using mainly the 12 CO emission in Section [ 3 ] 
The distances of these two cloud complexes are investi¬ 
gated in Section [4j Physical properties highlighting mass 
and mass surface density are discussed in Section [5| The 
dynamical nature of the clouds are investigated in Sec¬ 
tion [6] In Section [7| we discuss the method of predcting 
the SFRs and the evidence that the RCW 106 complex 
can be quanlified as a ministarburst. 

Finally, we summarize our findings in Section [8j 

2. DATA 

To investigate the properties of the region we used var¬ 
ious continuum and spectral line tracers, summarized in 
Tabled and discussed below. 


2 . 1 . 


12 CO molecular line from the ThrUMMS survey 

The ThrUMMS survey (jBarnes et al.||2013| covers the 
Galactic plane between 300° < l < 360" and — 1° < b < 
1° in the 12 CO (1-0), 13 CO (1-0), C 18 0 (1-0), and CN 
spectral lines at frequencies of 115.27, 110.20, 109.78, and 
113.60 GHz, respectively. [^] The survey was conducted 
using the new Very Fast Mapping (VFM) tech niqu e de¬ 


veloped for use with the 22 m Mopra telescope XV1 Half- 


Nyquist sampling rates were obtained along the scanning 
direction and between scan rows, yielding a 72" resolu¬ 
tion. Compared to the traditional on-the-fly mapping 
technique, the VFM technique is capable of observing 
a larger area in the same amount of time if Nyquist- 
sampled maps are not needed. As well as being more 
efficient, it provides a nearly uniform sensitivity per unit 
area. The data were converted from antenna tempera¬ 
ture T ant to main-beam temperature, T m b, by dividing 
T an t by the main-beam efficiencies 7^ m b of 0.42 for 12 CO 
(1—0), and of 0.43 for both 13 CO (1-0) and C 18 0 (1-0). 

2.2. 21 cm continuum and H i line emission from the 
SGPS survey 

From the Southern Galactic Plane Survey (SGPS; 


McClure-Griffiths et al. 2005) we extracted the H i 
atomic line" and 21 cm continuum data over the range of 
330° < l < 335° and b = ± 1°. The survey was conducted 
using the Australia Telescope Compact Array (ATCA) 

xv The data available to date for this study covered the range 
— 0?5 < b < 0?5. 

XV1 Operation of the Mopra radio telescope is made possible by 
funding from the National Astronomical Observatory of Japan, the 
University of New South Wales, the University of Adelaide, and the 
Commonwealth of Australia through CSIRO. 














































































Molecular cloud and star formation in the surrounding of the giant H n region RCW 106 
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Figure 1 . 21 cm continuum image (color scale) from the SGPS. Contour levels (white) are 0.2, 0.6, 1, and 2 Jy beam -1 . Within this 

region, we analysed continuum and molec ular li ne data listed in Table IT] in 31 distinct boxes (green). T he linear ruler ftop of diagram) 
assumes a distance of 3.6 kpc. OH masers \Caswell et al. 1980|) are plotted as red crosses, Jd II regions jjones & Dickey|2012) as magenta 
asterisks, and prominent OB clusters as blue stars. Further details are given in Appendix [Aland Table | 5 | 


interferometer and was complemented with the Parkes 
64 m telescope (FWHM = 15') for short spacings. The 
observations were performed simultaneously in a spectral 
line mode with 1024 channels across a 4 MHz bandwidth 
centered at 1420 MHz and in a continuum mode with 
32 channels across a 128 MHz bandwidth centered at 


Appendix [A] and Table [ 5 ] 

Plane Survey [MIPSGAL] | Carey et al.||2009 ). 


3. SPATIAL AND KINEMATIC STRUCTURE 
3.1. 21 cm continuum: defining the major structures 


1384 MHz (| McClure-Griffiths etahpOOlj 2005). For our 
study, we used data which have a FWHM 2.'2, line rms 
1.6 K kms -1 , and continuum rms < 1 mJy beam -1 (see 
Table [l]) . 

2.3. Infrared data from Herschel and Spitzer 

We used 160 and 250 /am images from the Herschel 
Space Observatory to examine more deeply the proper¬ 
ties of the cold dust near RCW 106. The fields were ob¬ 
served as part of the Herschel Infrared Galactic Plane 
Survey (Hi-GAL; Molinari et al. 2010|) at 70/160 /am 
with the P hotodetector Array C amera and Spectrome¬ 
ter (PACS;|Poglitsch et al.|2010|) and at 250/350/500 /im 
with the Spectral and Phot ometric Imaging REceiver 
(SPIRE; Griffin et al. 2010). The Hi-GAL data were 
taken in parallel mode with a fast scanning speed of 60" 
s -1 . The raw (level-0) data of each individual scan from 
both PACSand SPIRE were calibrated and deglitched 
using HIPE xvii | version 11.0. The SPIRE and PACS level- 


In the 21 cm continuum image (Fig. [T|) the giant H n 

region RCW 106 (Rodgers et al. 

I960]) is the largest 

(1?2 x 1°) and brightest (~ 287 J 
is an elongated structure in 332? t 

y) area. RCW 106 
)< l < 333?8 and 


—0?8 < b < 0?1, and angled roughly 60° with respect to 
Galactic aast within the /, b plane (hereafter referred to 
as the “aastern” part). In the region 330?5 < l < 331?7 
and —0?2 < b < 0?3 there is another similarly an¬ 
gled loose structure containing bright 21 cm continuum 
sources (hereafter the “western” part). Between these 
two complexes of thermal sources there is a gap around 
l rsj 332°. Projected in t his gap is an QB associa- 
tion R103B (332?08, -0?08) (iMel’Nik & EfremovfT995) 

A— o i ao a /ooiooa 


1 data were then fed to versi on 18 of the Scanamorphos 
software package xvm (Roussel 2013) to regrid and create 
the final maps. For our region, three Hi-Gal fields of 
2?5 x 2?5 were combined. The resultant 160 and 500 /am 
images have angular resolutions of 12" and 37", and la 
rms of 0.08 and 1.0 MJy sr -1 . 

To trace the warm dust emission associated with high- 
mass star formation, we used the 24 /im image from 
the Multi-band Imaging Photometer for Spitzer Galactic 

xvn HIPE is a joint development software by the Herschel Sci¬ 
ence Ground Segment Consortium, consisting of ESA, the NASA 
Herschel Science Centre, and the HIFI, PACS, and SPIRE consor¬ 
tia. 

xvm http://www2.iap.fr/users/roussel/herschel 


and also another further to the south, R103A (331V 
—0?93). These are at spectroscopic distances of 3.22 and 
3.00 kpc, respectively, and so in the foreground. The 
eastern part is located at a distance of ^ 3.6 kpc and the 
western part is located at a distance of ^ 5 kpc 
We have further divided the map in Fig. [l] into 31 sub- 
regions based on a contour corresponding to the 3cr noise 
level above the local background. Most of the sources are 
H ii regions whose properties are classified according to 
their morphologies in Appendix [Aj Additionally, there 
are four confirmed and two candidate supernova rem¬ 
nants (SNRs) ( [Green 2009). For full details of each sub- 
region refer to Appendix |A| and the summary provided 
in Table 02 

In Fig. [l] we also mark the positions of OH masers and 


H ii regions found in large s cale surveys (Caswell et al. 


1980| | Jones & Dickey |20 12 1). The positions of masers 
“Hmi 


and H n regions detected from radio recombination lines 
coincide with the bright radio continuum sources help 
to distinguish substructures and to quantify the cloud 
properties (see Sect.|4|. 


Jy beam 













































4 


H. Nguyen, Q. Nguyen Luong, P. G. Martin, et al. 



-140 -120 


-100 -80 -60 

Vlsr (km s- 1 ) 


40 -20 


Figure 2. The 12 CO spectrum averaged over the entire region 
(l = 330° — 335°, solid black), eastern part ( l = 332?5 — 334°, 
dotted green), and western part ( l = 330?5 — 332?5, dashed blue) 
as defined in Fig. [Tj 

3.2. CO and H i line emission 

To investigate the relative contributions of different 
clouds superimposed along the line of sight, we first plot 
the average 12 CO spectrum of the entire region along 
with the distinctive spectra of the two m ain s tructures, 
eastern and western, defined in Section |3.1| (Fig. [2]). 
In the average spectrum there are three distinct peaks 
within the velocity range —140 to 10 kms -1 . The spec¬ 
tra for both main structures exhibit significant peaks 
at roughly the same velocities, in the ranges [—112, 
—80], [—72, —60], and [—60, —35] kms -1 centered on 
Vlsr = —91, —63, and —48 kms -1 , respectively. These 
ranges are highlighted as green, blue, and red sectors in 
Fig-H For the western region, the three spectral max¬ 
ima are roughly the same, on average ~ 7 K. For the 
eastern region, there are two peaks in the lowest velocity 
range, the brighter one at Vlsr r\j — 55 kms 1 and the 
other at ~ —45 kms -1 (see Fig.j2|. The peak at ^ —55 
km s -1 coincides with the 13 CO peak that Bains et al. 


(2006) attributed to the cloud surrounding the RCW 106 
cluster. 

In Figs.[3|i-c and d-f we display maps of the 12 CO and 
H i emission integrated over three velocity ranges: [—112, 
—80], [—80, —60], and [—60, —40] kms -1 . Additionally 
in Fig. [3j the positions of OH masers and H n regions 
mostly coincide both spatially and spectrally with the 
molecular gas peaks. 

Fig. [3k (corresponding to the green velocity range in 
Fig. El) ahows extended CO emission spanning the entire 
longitude range. Its main peak is between 331° and 332° 
in longitude and extends just 0?2 below the mid-plane; it 
covers the western region and might extend to l = 333?4. 

The CO map integrated from —80 to —60 kms -1 (the 
middle white and blue velocity ranges in Fig. [2|, shown 
in Fig. ( 3 ) 3 , is dominated by the strongest emission from 
the Western region, but is otherwise relatively free from 
emission from the central plane (i.e., b ~ 0°). On the 
other hand, the map shown in Fig. H, integrated from 
—60 to —40 kms -1 (the middle white and red velocity 


ranges in Fig. D| is much brighter near the RCW 106 
H 11 region ananas a diffuse filamentary component in 
the western region. 

The H 1 emission in Figs. [3tl-f does not appear to cor¬ 
relate very well with either^ 2 COor 21 cm continuum 
emission, indicating that star formation activity does not 
correlate with this extended phase of the gas. However, 
in Figures [3]i-e, H 1 seems to form in the outer envelope 
of the MCCs. 

We have integrated the 12 CO, 13 CO, and H 1 data 
cubes in Galactic latitude (from —0?5 to 0?5) to pro¬ 
duce position-velocity (JL — V) diagrams (Fig. [4|. The 
main emission in the i — V maps is from —72 kms -1 to 
—35 kms -1 as in the averaged spectra (Fig.[2|. Another 
feature of strong though less prominent emission is from 
— 112 kms -1 to —80 kms -1 at low £ (western end). This 
is detached from the lower velocity structure supporting 
our division of the emission into two different complexes: 
in the East is the RCW 106 MCC spanning the velocity 
range from —72 to —40kms -1 , the West is the complex 
MCC 331.0+0.0 (V LSR = —90 kms -1 ) (hereafter referred 
to with the informative name MCC331-90 ) spanning the 
velocity range from —112 to —80 kms -1 . 

3.3. Insight from thermal dust emission 

Figure [5] shows images of the thermal dust emission 
obtained in the mid-IR by Svitze r an d in the far-IR and 
submillimeter by Herschel (Sect. 2.3) xlx These provide 
insight into the dust content, cloud morphology, and ra¬ 
diation field. In particular, data at 24 /im (Fig.[5k) allow 
us to trace local dust heating from the UV photons from 
young OB stars. Far-IR and submillimeter data probe 
cooler dust as well, both in the cloud and outside, and 
away from the influence of OB stars (Figures [Sb-c). 

The morphology of far-IR emission is generally similar 
to the CO emission (Fig. [3| suggesting that the gas and 
dust coexist within the MCCs. We also detect infrared 
dark clouds (IRDCs), seen as absorption in the mid-IR 
data but as emission in the far-IR data. These IRDCs 


are potential sites of massive star formation (Hennebelle 


et al. 20011; 

Rathborne et al. 

2006; bimon et al. 2006! 

Peretto et a 

1. 2010't Nguyen 

mong et al. 2011a). In 

the dust emission at all wavelengths we 
concentrations of dust at 331° and 333° 
morphology matches well with that in th 
tures seen in the radio continuum and i 
correlation indicates that both the 21 cm 
the dust heating are generated by the mas 
RCW 106 and MCC331-90 MCCs. The , 
data also reveal finger-like protrusions at 
the structure linked to MCC331-90. Th< 
plausibly formed by radiation pressure : 
OB stars in the local interstellar medium ( 

} see two main 
whose detailed 
e angled struc- 
u 12 CO. This 
continuum and 
sive stars in the 
Herschel far-IR 
; the bottom of 
sse features are 
from numerous 

[e.g., Krumholz 

& Matzner 2009; Critschneder et al. 2010; Tremblin et al. 


2013). 

In summary, the dust emission data from 8-500 /am 
confirm that this region has two distinct MCCs that are 
both extremely active in forming massive stars, including 
hosting a young stellar cluster. 

X1X Emission from Polycyclic Aromatic Hydrocarbons (PAHs) 
at 8 /im (not shown here) is also prominent in massive star forming 
regions (Peeters et al. |2004j ). 













































Molecular cloud and star formation in the surrounding of the giant H n region RCW 106 


5 


1 000 arcsec 


1 2CO + maser 4HII ^OB cluster 


-112 to -80 km s 





0.4 

- 0.2 

D 

J 

3 0 

J—0.2 

-0.4 


100 


100 


150 


100 

50 


0 

331 330 

— 1 1 2 to —80 km s 1 


334 


1 2C0 


333 Longitude (°) ^32 

4- maser *HII tyOB cluster 


— 80 to —60 km s 


-60 to -40 km s 


0 

200 


0.4 

0.2 

0 

0.2 

0.4 


0.4 

0.2 

0 

0.2 

0.4 


335 

1 000 arcsec 



-80 to -60 km s 



335 334 333 332 331 330 

Longitude ( ) 


Figure 3. Upper - 12 CO integrated map for (a): —112 to —80 kms 1 ; (b): —80 to —60 kms 1 ; (c): —60 to —40 kms 1 . The black 
contours are from the 21 cm continuum emission in Fig. [^with levels of 0.2, 0.6, 1, and 2 Jy beam -1 . As in Fig. ^positions of cataloged 
OH masers and H II regions are marked, but now only m the map with the appropriate velocity range. OB clusters are marked on all 
panels. Lower -Hi integrated maps ((d), (e), and (f)) for the same velocity ranges. The black contours are the 12 COemission from the 
corresponding panels above with levels in increments of 20 Kkms -1 from 10 to 90 Kkms -1 for (d) and (e) and 10 to 150 Kkms -1 for (f). 
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Figure 4. Hi (color), 12 CO (black contour, upper panel), and 13 CO (blue contour, lower panel) position-velocity diagrams 
integrated from b = — 0?5 to b = 0?5. As in previous figures, red crosses represent OH masers and magenta asterisks represent H II regions. 
The brightest 21 cm continuum sources in this region are at l = 331?5 and l = 333?6 (see Fig. [Hand Fig.pfi). H I absorption against these 
sources results in vertical bands of relatively low (green) net H I emission cutting through higher intensity (yellow, red) regions for the 
entire velocity range corresponding to foreground gas. 
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Figure 5. Dust emission (colour) from Spitzer 24/im (Top), Herschel 160 pan (Middle), and Herschel 500p,m (Bottom), overlaid with 
21 cm continuum emission contours (black). Blue stars mark the locations of bright OB clusters. 
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4. DISTANCES AND LOCATION OF THE OF RCW 106 AND 
MCC331-90 COMPLEXES 

As shown in Fig.[2j the RCW 106 molecular cloud com¬ 
plex has two peaks centered at —48 and —63 kms -1 , 
whereas the MCC331-90 molecular cloud complex has 
only one V LSR peak at ^ —90 kms -1 . We calculate the 


kinematic distances 

Reid et al. 

(2009), 

Duval et al. 

12009) 


allowing the approach of Roman- 
For RCW 106 and MCC331-90, 


respectively, in contrast to far distances 10.2-12 kpc and 
8.5-10 kpc, respectively. The far-distance solutions are 
very unlikely. First, at a distance of 11 kpc the MCC 
dimension would be ~ 200 pc,tripled the actual size cal¬ 
culated from the near solution. Second, the mass of 
RCW 106 would be a ~ 10 7 — 10 8 M 0 , as massive as 
the CMZ. Furthermore, at a Galactic latitude of -1° it 
would imply that these MCCs are ~ 300 pc off the Galac¬ 
tic mid-plane, which seems implausible even for the thick 
disc. The near kinematic distance to the RCW 106 com¬ 
plex is also consistent with the photometric distanc e es- 


timates to the RCW 106 OB cluster of 3.6 kpc (Moises 

eFalpOlT ). 


The 11 OH maser emission found by | Caswell et al.| 
(1980) he in the spatial and velocity range of either the 
ROW 106 or MCC331-90complex (see Figs. [II[3] and [4]), 
except for one source that lies in the spatial range of 
RCW 106 with an offset of ~ 40 kms -1 in V LSR . The 
average near kinematic distances of the OH masers also 
coincide with those of these two MCC, therefor e con - 


firmjDur argum ent for near kinematic distance. (Jones 
& Dickey |201 2) resolved the distance ambiguities for the 
RCW 106 and MCC331-90 complexes using H i absorp¬ 
tion and found that more than 75% have near distances. 
. Among these, 15 are assigned the near distances and 
four are uncertain. 

To be consistent with the measurements from other 
methods in t he literature (i.e., Bains et aL||2006 Moises 


et al 
RCW 


2011), we adopt 3.6 kpc as the distance to 
06. We also conclude that separation into two 


main molecular structures is physically meaningful, with 
both connected to signposts of massive star formation 
such as OH masers or H n regions. Hereafter, we adopt 
a distance of 3.6 kpc to the RCW 106 complex and 5 kpc 
to the MCC331-90 complex. 

With our calculated distances for these MCCs, we can 
infer their locations in the Milky Way (Fig. [6|. B ased 
on models of the Galaxy (George lin fc Ge orgelin 1976 


Rodriguez-Fernandez &; Combes [2005}" Reid et al.||2009p 

RCW 106 can be placed in the nearby Scutum-Centaur us 

arm, the major arm of the Galaxy, while MCC331-90is 
in the less prominent Norma arm. 

The Scutum-Centaurus arm, a counterpart to the 
Perseus arm (Fig. [6|) , is characterized by a high frac¬ 
tion of dense gas ([Sakamoto et al.||1997[ iRusseil et al. 


2005) which may be the reason why the most massive 


Young Massive Clus ters (YMCs) such as Westerlund 1, 
RSGC 1, RSGC 3 (|Portegie ~Z wart et ah 2010 l), and 


the ministarburst W43 (INguyen Luong et al. 2011b) are 


found therein. The £ — V diagram of the CO and H n 
emission shows a gradient toward the RCW 106 H n re¬ 
gion (see upper dashed white line in Fig. [I| that might 
reflect a velocity gradient along the Scutum-Centaurus 



Figure 6. Artist’s rendition of the Milky Way seen face-on by 
Robert Hurt of the Spitzer Science Centre with advice from Robert 
Benjamin from the University of Wisconsin-Whitewater. Colored 
dots mark the positions of GMCs, except for the black dot which 
marks the Galactic center. The RCW 106 complex is at l = 333? 0 
and has a distance of 3.6 kpc from the Sun, while MCC331-90is 
at l = 331? 8 and a distance of 5 kpc. Here it becomes obvious that 
the two GMCs are separate and actually part of different spiral 
arms. 


arm toward the RCW 106 complex. 

Although the Norma arm, a counterpart to the Sagit¬ 
tarius arm (Fig. [6]), has less active star formation overall, 
MCC331-90 is known to host the luminous and efficient 


massive star forming cloud G331.5-0.1 (Merello et al. 
2013a|b ). 


5. MASSES AND MASS SURFACE DENSITIES OF THE 
RCW 106 AND MCC331-90 COMPLEXES 

We obtain the mass from VF( 12 CO), the 
12 CO spectrum integrated over the velocity range —80 
to —40 kms -1 for RCW 106 and —112 to —80 kms -1 for 
MCC331-90. The molecular hydrogen column density 
is from 7Vh 2 = X x VP( 12 CO), with X = 1.8 x 10 20 
cm -2 K _1 kms -1 (Dame et al. 2001). This X factor is 


close to the value of 1.9 x 10 2U cm Z K 1 kms 1 derived 
for the Perseus molecular arm from the diffuse gamma- 


ray emission (Abdo et al. 2010) and is lower than the 
value of 2.75 x 10 2U cm _z K Tins -1 used for W43 and 


Cygnus X in 

Schneider et al. (2006) and Nguyen Luong 

et al. 

(2011b 

. Although A is uncertain because of the 
1 , metallicity, and excitation conditions, 

optical deptl 


the mass estimates are probably accurate to a factor 
of two. The mass M tota i in an area A c i ouc j would then 
be Mtotai = Wi 2 x Adoud x ju H2 m H where juh 2 = 2.8 is 
the mean H 2 molecular weight and mn is the H atomic 
mass. 

We use the 12 CO integrated intensity maps to define 
the extent of each cloud and adopt a rectangular shape 
that covers the main extent of the cloud (see Fig. pi) 
which is then used to calculate A using the assumed kine¬ 
matic distance. The areas of the two clouds are 2.6 x 10 4 
pc 2 and 2.2 x 10 4 pc 2 , which yield effective diameters of 
183 and 137 pc, respectively. However, we integrate only 
within the contour 7Vh 2 = 5 x 10 22 cm -2 (corresponding 
to W( 12 CO) = 277Kkms _1 ) to minimize the effects 
of the foreground and background emission. This yields 
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Table 2 

Properties of the two molecular cloud complexes in comparison with other star-forming complexes 


Complex 

-A c 1 0 U '.1 

(PC 2 ) 

D h 

(pc) 

Vel. range 0 
( km s -1 ) 

d d 

(kpc) 

M e 

(M 0 ) 

y f 

^gas 

(M 0 pc" 2 ) 

<T S 

( km s -1 ) 

M vir h 

(M 0 ) 

Q; vir 1 

^SFR j 

(M© yr -1 kpc -2 ) 

RCW 106 

2.6 xlO 4 

183 

-80 to -40 

3.6 

5.9T1.77 x 10 6 

221.2 

4.5 

2.1 x 10 6 

0.35 

0.15-0.3-(0.3-2.3) k 

MCC331-90 

2.2 xlO 4 

167 

-112 to -80 

5.0 

2.8T0.84 x 10 6 

127.1 

3.7 

1.4 x 10 6 

0.5 

- 

<Gould Belt> 1 

32.4 1 

6.4 1 

- 

0.27 1 

3.0 x 10 3 

79.3 1 

- 

- 



W43 m 

1.5X10 4 

~ 140 

80 to 100 

6.0 

7.1 x 10 6 

473 

9.3 

2.4x 10 7 

0.3-0.5™ 

0.65-(6.5) k 

Cygnus X m ’ n 

2.0X10 4 

~ 160 

-10 to 20 

1.7 

5.0xl0 6 

250 

4.2 

5.5xl0 6 

0.9 


CMZ m, ° 

9.6X10 4 

~ 350 

-225 to 225 

8.5 

3.0xl0 8 

3125 

- 

- 

- 



'Main velocity range of the structure. 
cm -2 K -1 kms — 1 except for Gould Belt 


“Surface areas calculated at adopted distances. ^Equivalent diameter D = a/ A c \ OVl ^/'k. 
d Adopted distance. “Mass calculated from the CO emission using an X of 1.8 X 10 20 

clouds being calculated from dust extinction map and for CMZ from thermal dust emission. -^Mass surface density E gas = M/A c i ou d- 
9 Velocity dispersion a = AVpwHM/v / 81n2. h Virial mass. *3D virial parameter described by ot v i T = M v i r /M tota i. J Star formation 
rate (SFR) density. k First numbers are the past SFR density and bold numbers in brackets are the future SFR den sity (see Sect.|7|. 
^Averaere of 20 lare-e molecular clouds from Spitzer cores to disks and Gould Belt surveys_jHejderman_et akj2 010[ ). m From [Nguyen] 
with X factor scaled to 1.8 X 10 20 cm -2 K -1 kms -1 . n From ISchneide^^^d*? I2006|) using 13 CO. “From 
using C ls O. 


Luone; et al. ( 

,2011b 

jUahmen et al.] 



masses 5.9 x 10 6 Mq and 2.8 x 10 6 M 0 and approximate 
surface densities 220 and 130 M e pc _2 for the RCW 106 
and MCC331-90 complexes, respectively. We estimate 
an uncertainty of 30% because of the background and 
foreground confusion as well as optical depth. The prop¬ 
erties are collected in Tabled 


of 13 CO (Bains et al. 2006), mm continuum ( 

Mooker- 

jea et al. 2004), and far-lR ( 

Karnik et al. 

2001 

toward 


sons. First, the previous measurements estimate only 
the clump mass whereas we calculate the total gas mass 
which includes in addition the more diffuse large-scale 
structure]^] Second, while these three studies focus in 
on the ^ 1 square degree area surrounding the RCW 106 
cl uster, we are studying a region almost five times larger. 

Heiderman et al. ([2010) investigated 20 low-mass star¬ 
forming clouds in the Gould Belt region and found an 
average mass^ 3000 M©, area 32 pc 2 , and surface den¬ 
sity of 8OM 0 pc -2 (Table §• Compared to these, the 
two complexes studied here are ~ 100 times larger in 
area and yet because the masses are so much higher the 
mass surface densities are also higher. In reference to a 
typical massive star forming region such as Cygnus X, the 
mass surface density of RCW 106 is comparable, while 
MCC331-90is about half. But compared to the mass 
surface densities of the extreme massive star forming re¬ 
gions such as W43 or CMZ molecular cloud, the values 
here are lower still (Table [2|. 

6. DYNAMICS OF RCW 106 AND MCC331-90 COMPLEXES 

We fit a single Gaussian profile to the spectra of 
the RCW 106 and MCC331-90 MCCs and obtains 
FWHM widths of 10.6 and 8.7 kms -1 corresponding 
to a one-dimensional (ID) velocity dispersion a = 
AWwhm/ V 81n2 of 4.5 and 3.7 kms -1 , respectively. 


in other observations (W43, 

Nguyen Luong et al. 

2011a 

Cygnus X, 
molecular cl 

Schneider et al. 

2006) or in simulations of 

iouds in spiral arms (Dobbs & Fringle 2013). 

First, to compare with other complexes, we estimate 


the gravitational instability criterium for the RCW 106 

xx The ratio of total clump mass to the total gas mass can 
also be considered as the clump formation efficiency, as elaborated 
further in Section [7] 


complex using the spherical approximation and assuming 
that its effective radius R sp h (half the diameter D in 
Table |2]) as: 


= 5 x 


-^sph^iD 

GM 




Mtotai 


(i) 


Both complexes have low spherical virial parameters ( 
0.35 for RCW 106 and 0.5 for MCC331-90 ; see Table [2|. 

However, since the global structure of RCW 106 com¬ 
plex is filamentary with an aspect ratio r = A = = 


3.6, we use the gravitational instabilit y criter ium of an 
infinite-length isothermal gas cylinder (Ostriker 1964) In- 
utsuka fe Miyama||1992| see also IDIBaffl 1958710 zernoi 
1964p . This model defines a critical line mass 


mass 


per unit length, Mii nejCr i t solely on the sound speed c s , 
above which the filamentary cloud is unstable against 
collapsing. For a large molecular cloud complex, the in¬ 
ternal motion depends also on turbulence in addition to 
its thermal motion, therefore we assume that the appar¬ 
ent ” sound sp eed c,% equals th e observed line width AV 
(similar as in Dobashi et al. |2014 ). We calculate the 

-f-l 7 n nm o nn -h n i nnnnr'nnm oionr* /'A/Pl r\ va od o o • 


filamentary virial parameter a ; 
Mine ,crit _ 2 cl/G 


a m = 


defined as: 
465AR [Mqpc- 1 ] 


^line, 


gas 


M/l M/l [Mqpc- 1 ] 


( 2 ) 


The total gas mass M is integrated over the RCW 106 
and MCC331-90 filamentary clouds with length l and 
width w. We obtain low <a 0r of 0.15 for RCW 106 and 
0.26 for MCC331-90 , which confirm that these two com¬ 
plexes are indeed gravitationally unstable. 

While the gravitationally unstable states of these 
clouds are certain, a more detailed investigation is needed 
to justify their origin. They might be caused by strong 
dynamical affects related to their locations spread along 
the respective spiral arms (see Section [ 4 ]) and/or feed¬ 
back from the vigorous burst of star formation within 
(see Section [7|). 

7. A BURST OF STAR FORMATION IN THE RCW 106 
COMPLEX 

Lyct continuum photons from young stars create H 11 
regions quantifiable by free-free emission. We find 30 
H 11 regions hosting O- or B-type stars in the RCW 106 
complex (Table [ 3 ]). We checked the CO and H 1 R LS r of 
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gas surrounding these H n regions and confirmed that 
they are in the velocity range of the RCW 106 complex. 
We estimated the radio continuum flux density S v of each 
H ii region by aperture photometry on the the 21 cm map 
within the contour defined by the 3cr noise level (see Ap¬ 
pendix 0 . An average background level of ^ 0.09 Jy 
beam -1 was derived by averaging the intensity wit hin 
a surrounding empty region. Following Mezger & Hen¬ 
derson (1967), the number of Lya continuum photons 


powering an ionization-bounded region can be computed 
using 


Nhya 

_ & 

V 

0.1 

T e 

-0.45 

" d " 

8.9 x 10 46 s -1 

Jy 

GHz 


10 4 K 


kpc 


,__ (3) 

where T e = 8000 K (|Wilson et al. 2012) is the electron 


temperature, v — 1.42 GHz is the observing frequency, 



spectral type 07 or earlier in each H n region of the 
RCW 106 cloud complex (Table [3]); in total, there are up 
to 54 07 stars currently formed. Assuming that the stel¬ 
lar masses in the RCW 106 complex follow the Salpeter 
Initial M ass F unction (IMF) dN/dm = A x m -2,35 


Salpeter 1955) and have a minimum cut-off mass of 
.08 M© and maximum cut-off mass of 50 M 0 , the to¬ 
tal stellar mass of RCW 106 will then be 


AL 


-d 


—0.35 _ m -0-35 

m- a dm = Ax —*3“-siin_ 

—0.35 


( 4 ) 


The normalization factor A can be calculated from the 
total number of 07 stars calculated above by the approx¬ 
imation that this is the total number of all stars with 
masses in the range of 25 to 50 M e so that 


N" 


A x 


rn 


-1.35 


— rn 


-1.35 


-1.35 


( 5 ) 


To gauge the uncertainty, we varied the upper mass 
limit to 100 M 0 and allowed for a dispersion of ±5 stars in 
the estimate of the total number of 07 stars. Finally, we 
obtain a value of ~ 776ll^g 3 for the factor A and a total 
stellar mass of 48q 7 x 10 3 M 0 . Replacing the Salpeter 
IMF by the Kroupa IMF (|Kroupa 200l]) lowers the total 
stellar mass but it is within the calculated uncertainty. 

7.1. The immediate past 

Because an 07V star with a mass of 25 M 0 has a life¬ 
time of - 2 x 10 5 yr starting from its first appearance 
on the main-sequence diagram, we can calculate the 
star formation efficiency e£^ d , the star formation rate 

XX1 Of course, it is possible that the ionization is produced by 
a combination of several stars of somewhat later spectral type; 
however, usually only the brightest star dominates the output of 
ionizing photons. 


SFR^ and the SFR density cloud of the entire 
RCW 106 cloud complex for this “past” time interval as 


past 

^cloud 


Af* 


Afcloud ± ^ 


= 0.008t°;°03, 


o p DPast _ Af* 

cloud 2 x 10 5 yr 


— 0.25^q 023 M 0 yr 


-l 


y^past 

^SFR, cloud 


SFR 


icloud 


= 9.5to;g M© yr 1 kpc 


-2 


( 6 ) 


( 7 ) 


( 8 ) 


with the total stellar mass M * = 4.8 x 10 4 M 0 , the 
total gas mass M c \ ou d = 5.9 x 10 6 M 0 calculated over an 
area of A c \ ou a = 2.6 x 10 4 pc 2 . 

The estimated e^ s ^ d is in the middle of the range for 
GMCs, 0 . 002 - 0 . 2 , as deri ved from the ionizing flux of 
young stars ([Murray] 2011). However, the SFR den¬ 
sity of RCW 106 is high and for its surface density of 
220 M 0 pc -2 (Sect. p| it is well above the trend in the 
Schmidt-Kennicutt relation (see Fig. [7]). This shows that 
RCW 106 has been very active in forming massive stars 
during the last 2 x 10 5 yr. 

7.2. The immediate future 

The large amount of gas remaining in this complex 
suggests that the burst of star formation activity in the 
RCW 106 molecular cloud complex might not yet be fin¬ 
ished. We investigated this by estimating the future star- 
formation activity using the dense clump and core pop¬ 
ulations. 

First, we make an assumption that the dense clump 
mass of 2.2 — 2.7 x 1 O 5 M 0 derived from 13 CO line emis- 
sion (|Bains et al7|| 20Q6 ) or 1.2 mm continuum emission 
(Wong et al. 2008} represents the clump mass of the en¬ 
tire RCW 106 "complex. This assumption is justified since 
most of the dense gas conc entrate s in the region co nsid- 
ered by Bains et al.|j2006) and Wong et al. (2008]) (see 
for example Figure [5]) . W e est imate the clump forma¬ 
tion efficiency, c c i 0 udiciump 


as the ratio of the total 


clump mass to the total cloud mass (Eden et al. 2012 


2013 


Lo uvet et al.||2014 ): 




^cloud—>-clump — 


clump 


A7 c loud 


0.042 ±0.013. 


( 9 ) 


The estimated e c i ou diciump of the RCW 106 complex is 
about 10 times higher than that of the famous massive 
star-forming regio n Cygnus X an d almost equal to that 
of W43 (Nguyen Luong et al. 2011b). Thus the RCW 106 
complex is forming massive clumps almost as efficiently 
as the W43 star forming region despite the mass surface 
density of RCW 106 being four times lower. 

Secondly, a cloud’s efficiency at converting mass into 
dense clu mps subsequently affects the star form ation ef¬ 
ficiency (Eden et al. 2012 ; Louvet et al.| 2014|). If the 
dense clumps [r ~ 1 pc) in the RCW 106 complex form 
dense cores (r ^ 0.1 pc) with a mass transfer efficiency 
^clumpi-core 5 and dense cores will form protostars with a 
mass transfer efficiency e cor ei*? the future star formation 
efficiency of the entire cloud will be 

^cloudi* = ^cloudiclump X ^clump—>-core ^ ^corei* • (10) 

xxn This paramet er is also known as t he compactness of a molec¬ 
ular cloud ([Nguyen Luong et al. 2011b). 
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Figure 7. The Schmidt-Ken nicutt relat ion between the star formation rate (SFR) density and gas surface density extending from normal 
spiral galaxies to starbursts (|Kennicutt|[l998|). Also plotted are values for the average Milky Way, the immediate past and future of the 
RCW 106 complex, the past and future SI R density of the W43 complex, the past and future SFR density of the Cygnus X complex, and 
the past SFR density of the CMZ complex. 

galactic plane within a few degrees of the Galactic centre, 
well-known for its inefficiency in star formation despite 
being 100 times more massive than RCW 106 complex 
Kruijssen et al.|2014 ). 


(Parmentier & Pfalzner 

2013 Louvet et al. 2014) and 

^core—>•* 2L- 

0.3 - 0.¥ ( 

Alves et al. 2007 Parmentier 

& Pfalzner 

2013 Myers 2014|), we will have a future 


(Immer et al. 2012 


of 


Because the m assive protostellar phase has a timescal e 
2 x 10 5 yr ( Motte et al.||2007 Russeil et al.||2012 ), 


(“future”) SFR and SFR density as 


yr 


SFR 


future 

cloud 


^cloud—>•* 

2 x 10 5 yr 


Mcloud = 0.12 ± 0.1 M e yr - 


and 


(ii) 


future 

^SFR, cloud 


C rp zofuture 

- = 4.8 ± 3.8 M 0 yr _1 kpc -2 , 


While this is consistent with the large mass reservoir 
still available in the RCW 106 complex, there are still 
questions outstanding including what the final conver¬ 
sion efficiency will be and its relationship to the virial 
parameter (Section [6|. 

7.3. Uncertainties 

The calculations in Sections rm and 17721 should involve 
the general uncertainties from the assumptions of the 
IMF, the timescales, and the mass transfer efficiencies. 
Several authors suggest that the IMF of the extr eme star 


4 


cloud 


formation environment is t op-heavy (Stolte et al. 


( 12 ) 

This SFR density and the surface density also combine 
to place the RCW 106 complex in t he starbu r st reg ime 
of the_ Schmidt-Kennicutt relation (Kennicutt|[l998; see 
Fig.0. It should be possible to refine and confirm the 
future SFR density of the RCW 106 complex by count¬ 
ing the massive dense cores identified using the Herschel 
data. This will be the topic of a subsequent paper. 

These SFR densities are similar to the future SFR den¬ 
sity of W43 and higher than that of Cygnus X, reinforcing 
the conclusion that RCW 106 is undergoing a ministar- 
burst at present. The similarity between past and future 
SFR densities implies that this ministarburst event is not 
yet finished but is in the declining phase. It is also much 
higher than the SFR density of CMZ, the region of the 


2005 Chabrier et al. 

2014J) 

this statement is still 

Largeb 


Howeve r, the valid ity _o. 


2002 


2006 


eeper 

slope, the SFR likely increases by about 20%. We as¬ 
sume a burst period timescale of 2 x 10 5 yr, which might 
have variations within a factor of two; but this is the 
most appropriate assumption since this is the measured 
timesc ale of an OB star (Martins et al. 2005) and an 
MDC (|Motte et al.||2007|). The third uncertainty is the 
mass transfer efficiencies from clouds to clumps, to cores 
and to stars. Our assumptions of these quantities are the 
average results of the current state-of-the-art theoretical 
and observational studies, therefore the true values might 
vary at most 10%-30%. We already take these three un¬ 
certainties into account in our calculations. 

The uncertainties from measurements of Ly<a photons 
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and the clump mass are low since we focus only on the 
brightest regions in the radio continuum map (in case 
of estimating the Ly<a photons) and in millimeter con¬ 
tinuum or 13 CO maps (in case of estimating the clump 
mass). For the later case, the clump mass is verified with 
an in depen dent measurement from NH3 emission (Lowe 
e t~aLl|20l4 ). 


8. CONCLUSIONS 


Using a combination of data sets, we have defined and 
characterized two MCCs in the general direction of the 
RCW 106 OB cluster, and we have labelled them the 
RCW 106 and MCC331-90 complexes. Their properties 
are summarized as follows 


• The U LSR for gas in the RCW 106 complex ranges 
between —40 and —80 kms -1 and the U LSR for 
MCC331-90 is between —80 and —112 kms -1 . Both 
complexes have velocity dispersions comparable to that 
of Cygnus X (cr=4.2 kms -1 ) but smaller in comparison 
to W43 (<t= 9.3 kms -1 ), which is right in the Galactic 
bar and a known turbulent area. 

• OH masers and H 11 region studies confirmed the near 
kinematic distances of 3.6 kpc for RCW 106 and 5 kpc for 
MCC331-90. This places the RCW 106 molecular cloud 
complex in the Scutum-Centaurus arm and MCC331-90 
in the Norma arm. 

• The RCW 106 complex has a mass of 5.9 x 10 6 M© 
and a surface density of ~ 220 M© pc -2 , whereas the 
MCC331-90 complex has a mass of 2.8 x 10 6 M© and 
surface density of ^ 130 M© pc -2 . These surface den¬ 
sities are higher compared to the average Gould Belt 
cloud (^ 70 M© pc -2 ) and the surface areas are 10 3 times 
larger. 

• These two complexes are of comparable size 
(d ~180 pc) but have less mass in comparison to other 
large GMCs such as W43 and Cygnus X. Hence, the sur¬ 
face densities are smaller. 

• The virial parameters are greater than unity, indicat¬ 
ing that they are gravitationally unbound. 

• Using the 21 cm continuum, we separated the region 
into 31 subregions: 25 containing H 11 regions and 6 con¬ 
taining SNRs. We estimated that there are about 50 
young 07V stars currently in the RCW 106 complex. 

• For the RCW 106 molecular cloud complex, we derive 
a past global star formation efficiency, SFR, and SFR 
density and estimate values to quantify star formation in 
the near future. These values suggest that the RCW 106 
complex is undergoing a ministarburst event. 
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APPENDIX 

A. 21 CM CONTINUUM SOURCES 

In Fig. [l] we selected 31 rectangular subregions in the 21 cm continuum data in which there is significant emission. 
These are the “boxes” whose details are captured in Table [3] 

Across the continuum image the background is rather similar, about ^ 0.09 Jy beam -1 , and contours are given 
starting at the 3 a noise level above the local background. We integrated the background-subtracted 21 cm intensity 
within the area defined by the 3cr noise contour. The flux densities range from 0.1 to 12 Jy and the sizes range from a 
few arcminutes to several hundred arcminutes. There is no apparent correlation of size and flux density. From these 
parameters and the distance, we have derived and tabulated the characteristics of the star formation in Table [3j Star 
formation in the RCW 106 complex is discussed in Section [7| 

Because of the massive star formation across the region there are many H n regions. These can be classified according 
to their shapes: circular or irregular. The H n region around RCW 106 is classified separately as a giant H n region 
and the components are displayed in Figure [8j There are also several SNRs seen in projection on the region. These 
structures can be seen in magnified views of these “boxes” in Figs. [9}|TT]in the subsections below. Within each figure, 
the boxes are shown on the same angular scale, as labeled, with the box number at the lower left of each sub-image. 
The spatial locations and relationships of the boxes can be seen in Fig. [l] In all figures, except Fig. [Tljfor SNRs^a 
linear scale is also given, assuming a distance of 3.6 kpc, although the structures associated with MCC331-90 (Table[3| 
are at 5 kpc. 
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Table 3 

Details of boxes and subregions as selected from the 21 cm continuum map 


No 

i b 

n n 

Type 

Vlsr 

(km s -1 ) 

Cloud 

Ang. size 

O 

Area S h 
(pc 2 ) 

,(21 cm) 

(Jy) 

SFR density 
(M 0 yr” 1 kpc” 2 ) 

Mass 

(M 0 ) 

^gas 

(M© pc -2 ) 

log JV Lyo , 
(s- 1 ) 

# 07V 

M tot 07V 
(M 0 ) 

1 

334.7 0 

CH II 

-60 to -30 

RCW 106 

2.6 

- 

0.48 

0.73 

- 

_ 

47.68 

0 

0 

2 

334.7 - 0.6 

CH II 

-60 to -30 

RCW 106 

2.9 

- 

0.96 

1.3 

- 

- 

48 

0 

0 

3 

334.6 - 0.1 

CH II 

-60 to -30 

RCW 106 

5.5 

- 

1.8 

0.78 

- 

- 

48.31 

0 

0 

4 

334.6 0.4 

CH II 

-60 to -30 

RCW 106 

1.6 

- 

0.12 

0.53 

- 

- 

47.21 

0 

0 

5 

334.5 0.8 

CH II 

-60 to -30 

RCW 106 

3.7 

- 

0.29 

0.14 

- 

- 

47.76 

0 

0 

6 

334.2 0.1 

SNR 

- 

- 

27.7 

- 

6.4 

0.24 

- 

- 

48.88 

- 

- 

7 

333.6 - 0.2 

GH II 

-97 to -78 

MCC331-90 

31.5 

2109.35 

- 

- 

1.7 xlO 5 

80.95 

- 

- 

- 




-58 to -32 

RCW 106 


881.48 

85 

2 

2.0x10 s 

227.99 

50.02 

20 

500 

8 

333.7 - 0.4 

CH II 

-80 to -57 

RCW 106 

2.1 

36.04 

0.23 

0.61 

2.7x 10 3 

74.15 

47.29 

0 

0 

9 

333.3 - 0.4 

H II 

-93 to -77 

MCC331-90 

16.6 

496.31 

38 

2.2 

2.2x 10 4 

44.45 

49.65 

- 

- 

10 

333 - 0.4 

H II 

-76 to -36 

RCW 106 

19.6 

402.26 

73 

2.6 

1.2x 10 4 

28.77 

49.94 

17 

425 

11 

333 - 0.7 

H II 

-76 to -36 

RCW 106 

12.9 

- 

4.9 

0.57 

- 

- 

48.74 

1 

25 

12 

332.9 0.8 

CH II 

-80 to -57 

RCW 106 

5.6 

- 

6.6 

2.1 

- 

- 

48.93 

1 

25 

13 

332.7 -0.6 

H II 

-60 to -30 

RCW 106 

16.8 

103.62 

36 

1.6 

l.OxlO 4 

99.17 

49.65 

8 

200 

14 

332.5 -0.1 

CH II 

-106 to -76 

MCC331-90 

4.3 

83.84 

1.9 

1.1 

9.8x 10 3 

117.17 

- 

- 

- 




-72 to -32 

RCW 106 


49.02 



3.6x10 s 

73.84 

48.37 

0 

0 

15 

332.4 0.5 

SNR 

- 

- 

14.6 

644.09 

18 

0.96 

5.8xl0 4 

90.73 

49.4 

- 

- 

16 

332.4 -0.4 

SNR 

- 

- 

9.7 

- 

20 

1.9 

- 

- 

49.41 

- 

- 

17 

332.2 -0.1 

CH II 

-101 to -77 

MCC331-90 

1.8 

49.64 

0.12 

0.53 

6.4xl0 3 

129.6 

- 

- 

- 




-77 to -48 

RCW 106 


32.21 



1.6x10 s 

50.44 

47.2 

0 

0 

18fl 

333.1 0 

CH II 

-98 to -82 

MCC331-90 

3.9 

748.75 

1.6 

1.4 

6.0X10 4 

80.66 

48.26 

- 

- 




-73 to -36 

RCW 106 


469.51 



3.8X10 4 

81.75 

48.26 

0 

0 

18f2 

- 

CH II 

-60 to -30 

RCW 106 

7.9 

- 

1.9 

0.49 

- 

- 

48.98 

1 

25 

18f3 

- 

H II 

-60 to -30 

RCW 106 

10.9 

- 

7.3 

1.01 

- 

- 

48.39 

0 

0 

19 

332.1 -0.4 

CH II 

-61 to -40 

RCW 106 

7.2 

65.68 



2.5x10 s 

38.61 

- 

- 

- 

20 

331.9 0.2 

SNR 

- 

- 

11.1 

368.16 

3.7 

0.56 

3.3X10 4 

89.6 

48.78 

- 

- 

21 

331.6 0 

SNR 

- 

- 

13.6 

1097.25 

41 

2.4 

1.8x10 s 

166.59 

49.71 

- 

- 

22 

331.4 -0.3 

H II 

-107 to -88 

MCC331-90 

9.7 

407.17 

7.9 

1.5 

4.3x 10 4 

105.16 

- 

- 

- 




-74 to -36 

RCW 106 

216.67 




3.1 x 10 4 

144.96 

48.97 

1 

25 

23 

331.3 0.5 

CH II 

-61 to -40 

RCW 106 

8.2 

- 

2.4 

0.59 

- 

- 

48.59 

0 

0 

24 

331.3 0 

CH II 

-104 to -72 

MCC331-90 

7.5 

154.31 

6.3 

1.3 

2.6xl0 4 

169.12 

- 

- 

- 




-72 to -37 

RCW 106 


96.66 



8.2x10 s 

84.42 

48.91 

1 

25 

25 

331.2 -0.2 

CH II 

-104 to -79 

MCC331-90 

4.9 

123.88 

3.4 

1.7 

2.3X10 4 

183.98 

- 

- 

- 




-68 to -32 

RCW 106 

69.47 




2.2x10 s 

32.31 

48.63 

0 

0 

26 

331.1 -0.5 

CH II 

-68 to -32 

RCW 106 

3.6 

- 

0.98 

0.69 

- 

- 

48.01 

0 

0 

26B 

- 

CH II 

-68 to -32 

RCW 106 

3.7 

- 

1.1 

0.84 

- 

- 

48.14 

0 

0 

26C 

- 

CH II 

-68 to -32 

RCW 106 

5.5 

- 

3.4 

1.3 

- 

- 

48.64 

0 

0 

27 

331.1 -0.2 

H II 

-101 to -76 

MCC331-90 

12.3 

481.05 

9.3 

0.95 

8.3xl0 4 

171.78 

- 

- 

- 




-69 to -33 

RCW 106 

286.99 




1.6xl0 4 

54.77 

49.09 

2 

50 

28A 

330.7 -0.4 

H II 

-101 to -80 

MCC331-90 

11.4 

528.18 

11 

1.3 

2.9xl0 4 

55.97 

49.15 

- 

- 




-69 to -32 

RCW 106 

332.42 




4.1 x 10 4 

122.17 

49.15 

2 

50 

28B 

330.7 -0.4 

CH II 

-69 to -32 

RCW 106 

3.9 

- 

1.8 

1.4 

- 

- 

47.54 

0 

0 

28C 

330.7 -0.4 

CH II 

-69 to -32 

RCW 106 

2.9 

- 

0.27 

0.56 

- 

- 

48.37 

0 

0 

29 

330.3 -0.4 

CH II 

-80 to -30 

RCW 106 

2.8 

13.68 



1.2x10 s 

88.77 

48.37 

0 

0 

30 

330.2 1 

SNR 

- 

- 

5.9 

244 

1.3 

- 

- 

- 

- 

- 

- 

31 

330.1 -0.1 

CH II 

-112 to -80 

MCC331-90 

3.2 

45.74 

- 

- 

5.3x10 s 

115.8 

- 

- 

- 




-67 to -33 

RCW 106 


23.44 

0.67 

0.86 

1.0x10 s 

43.74 

47.89 

0 

0 


A.l. RCW 106: a giant H n Region structure 

The H ii regi on complex RCW 106 spans approximately 1?5 in the range 332° < l < 333?5 and so at a distance of 
3.6 kpc (Lockman[l979) is of “giant” extent, 94 pc. The structure lies predominately below the Galactic mid-plane. 
We separated the structure into six components which are large themselves, illustrated in Fig. [8] and described in the 
brief notes below. 

Box 7 (RCW106a). Large in size at nearly 33 pc. It has many known H n regions, star clusters, and SNR G333.6-00.2 
which is located at the brightest peak (see Fig. [8k). 

Box 9 (RCW106b). An extended structure of ^ 17 pc in length (see Fig. It has an irregular shape with 

prominent protrusions. The main intensity peak is crescent shaped. 

Box 10 (RCW106c). An irregular H n region with a fairly large main peak and size of almost 9 pc (see Fig. [8J3). 
Box 11 (RCW106d). An irregular H 11 region with a fairly large main peak and a tail (see Fig. J8n). 

Box 13 (RCW106e). An irregular structure with three main peaks loosely connected and a size of about 17.5 pc 

(Fig- [§£)• _ 

Box 18 (RCW106f). Diffuse H 11 region with four dominant peaks (Fig.pt). 


A.2. Circularly shaped H 11 regions 

A number of the smaller H 11 regions in Fig. 1 appear to be fairly “circular” in structure. These are illustrated in 
Fig- E and described below. 

Box 1. The compact H 11 region contains the single IRAS source 16226-4900 at its center. 
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Figure 9. Circularly shaped H II regions. The distance scales use both near and far kinematic distances. 

Box 2. An ultracompact H n region coincides with an IR star cluster [DBS 2003] 170 


et al. 


(2003), the composition is a bit ambiguous but there are a few bright _stars 


As discussed by 
The Lt ii region GAL334.71 


Dutra 

TKWT 

found in tms box is 17.4 kpc away as determined by its positive velocity (Russeil et al.||2005|). 

Box 3. The compact H n region that we detect is rough ly 5.8 pc in size and probably coincides with the H n region 
GAL 334.68-00.11. It has a nearby cluster of IR sources (Ro bitaille et al.||2008 ). 

Box 4. This compact H n region is roughly 1.6 pc in size, which almost makes it an ultracompact candid ate. 

Box 5. There is a YSO on the border of the compact H n region. The star cluster [MCM 2005b] 79 (Mercer et al. 
2005) is also nearby. This source has a faint protrusion in the lowest contour, not so significant that we describe" it"as 


irregular. 

Box 8. This compact H n region has a few YSOs nearby. 

Box 12. This compact H n region contains the IR star cluster [DBS 2003] 102. It also coin cides with the known H n 


regions GAL 332.98+00.79 and [WHR97] 16112-4943. The latter is classified as UCHII (Walsh et al 
Box 14. This com pact g rouping encompasses the known H n regions [KC 97c] G332.5-0U 


i 


1997). 

54-00.11, and 


. . w L j I7GRRR 

GAL 332.54-00.11 (Kuchar & Clark] 1997) and the star clusters |DBS 2003] 160 and 161. This clump is a potential 
OB cluster as seen from the numerous OB stars within the region. There are also some YSOs and outflow regions 
within the clump. 


xxm [ref.]# from SIMBAD: http://simbad.u-strasbg.fr 
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Figure 10. Irregularly shaped H II regions. The distance scales use both near and far kinematic distances. 

Box 17. This is the compact H n region IRAS 16119-5048 driving an outflow. 

Box 19. An H n region hosting a star cluster [MCM 2005b] 74. It has a cometary shape and is significantly brighter 
with an off-center peak in its circular structure. It is located near the supernova remnant, SNR G332.4-0.4. 

Box 23. Most likely the compact H n region GAL 331.36+00.51, the structure is elliptical with a slight protrusion 
along the longitudinal axis and is 8.6 pc in size. 

Box 25. A possible compact H n region. It is ~ 5 pc in size with the known H n regions [KC 97c] G331.3-00.2, GAL 
331.28-00.19, and GAL 331.26-00.19 within its peak. 

Box 29. A compact H n region, 2.9 pc in size and associated with the IR cluster [DBS 2003] 153. It also hosts the 
know H ii region IRAS 16037-5223 and GAL 330.31-00.39. 

Box 31. A compact H n region probably powered by GAL 330.04-00.05. 


A.3. Irregularly shaped H n regions 

The remaining H n regions as seen in Fig. [10] are “irregular” in outline as compared those in Fig. [9] Descriptions of 
these follow. 

Box 22. An extended H n region of ~ 10 pc in size. It has numerous IR star clusters and known H n regions. It 
displays a comet-like structure with its brightest and largest feature at 331?33. There are two smaller peaks trailing 
off from the tail. They are diffuse and seemingly attached to the main structure. 

Box 24. A compact H n region with the IR star cluster, [DBS 2003] 158. It is elongated and has double peaks. 
The main peak is off-center in one of the protrusions where the known H n regions are located. It is also near the 
supernova candidate MSC331.8+0.0 whose structure may be related. 

Box 26. The structure is elongated with three main compact peaks that may all be connected. Together they form 
an elongated structure that is ^ 11 pc in length. 

Box 27. Both GAL 331.03-00.15 and PMN J1610-5150 coincides with this H n region. There are numerous OB 
associations spread around the main intensity peak. 

Box 28. This has been subdivided into three regions 28A, 28B, and 28C from largest to smallest with sizes ~ 12, 4, 
and 3 pc, respectively. The H n region 28A hosts numerous known H n regions and the IR star cluster [DBS 2003] 155 
while 28B appears as a compact H n region with IR star cluster [DBS 2003] 154. 28C is as yet undefined. Individually, 
28A has a comet-like shape with two tails while 28B and 28C are circular in nature. 


A.4. Supernova remnants and candidates 

There are four confirmed and two candidate supernova remnants (SNRs) (Fig. pT]). These are seen in projection and 
are not necessarily associated with the two main star-forming complexes. 

Box 6. The SNR Candidate G333.9+00.0 seems very faint with many separated clumps. 

Box 15. The SNR G332.4+0.1 has three main peaks that may be part of one shell judging from the contour lines. 
From OH absorption it is thought to be in the range of 7-11 kpc (Caswell & Haynes| [l975). 

Box 16. The SNR G332.4-0.4 seems to have a circular shell. A higher intensity portion situ ated near the bottom is 


2009 


Hi 


1 & Haynes|1975 ). 


the RCW 103 star cluster. It has a kinematic distance of 3.1 kpc from H i absorption lines (Green 
Box 20. The shell of SNR G332.0+0.2 is a crescent and is thought to be at at least 7 kpc (Caswell 
Box 21. The SNR candidate MSC 331.8+0.0 shows a h igh i ntensity portio n with scattered filaments? Box 3C 
SNR G330.2+01.0 has a minimum distance of 4 kpc (McClure-Griffiths et al.pO Ol). 


The 
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